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Study of t/>( 3770) decaying to Baryon anti-Baryon Pairs 
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Department of Physics, Tsinghua University, Beijing 100084, People’s Republic of China 

To study the decays of 0(3770) going to baryon anti-baryon pairs ( BB ), all available experiments 
of measuring the cross sections of e + e _ —> BB at center-of-mass energy ranging from 3.0 GeV to 
3.9 GeV are combined. To relate the baryon octets, a model based on the SU(3) flavor symmetry 
is used and the SU(3) breaking effects are also considered. Assuming the elctric and magnetic form 
factors are equal (|G_b| = |Gm|), a global fit including the interference between the QED process 
and the resonant process is performed. The branching fraction of ip( 3770) —> BB is determined to 
be (2.4±0.8±0.3) x 10~ 5 , (1.7±0.6±0.1) x 10“ 5 , (4.5±0.9±0.1) x 10“ 5 , (4.5±0.9±0.1) x 10~ 5 , 
(2.0±0.7±0.1) x 10 _s , and (2.0±0.7±0.1) x 10 -5 for B = p, A, £ + , £°, and E°, respectively, where 
the first uncertainty is from the global fit and the second uncertainty is the systematic uncertainty 
due to the assumption \Ge\ = |Gm|- They are at least one order of magnitude larger than a simple 
scaling of the branching fraction of J/ip/ip(3686) BB. 

PACS numbers: 13.25.Gv, 11.30.Hv 


I. INTRODUCTION 


The ip( 3770) is the lowest lying 1 charmonium state 
above the charmed meson pair threshold. It decays dom¬ 
inantly into D°D°/D + D~ while the decays to the light 
hadron (LH) final states are OZI-suppressed. It is still 
unclear about the nature of ip( 3770). If it is a pure cc 
bound state, the branching fraction of '0(3770) into non- 
DD decays ranges from less than 1% from the potential 
models mu to about 5% from the non-relativistic QCD 
calculations mm. If 0(3770) has a four-quark admix¬ 
ture, the total non -DD branching fraction could be up 
to 10% [5|. 

Experimentally, the BES collaboration reported a large 
non-HD branching fractions of (14.5 ± 1.7 ± 5.8)% [BH8] 
neglecting the interference between the 0(3770) resonant 
amplitude and the QED continuum amplitude. Only 
considering the interference between the one-photon am¬ 
plitude of the 0(3770) resonance and the QED continuum 
amplitude, the CLEO collaboration found this branching 
fraction to be (—3.3 ± 1.4l|;|)% |]J]. To clarify the dis¬ 
agreement, many exclusive non-DU decays with the light 
hadron final state have been searched for using two meth¬ 
ods PME3- One method is to compare the cross section 
at the center-of-mass (c.m.) energy (y/s) close to the 
0(3770) nominal mass and that far from any charmonium 
resonance (for example, the two energies are 3.773 GeV 
and 3.671 GeV for the CLEO collaboration). Only for 
the final state (pr], there is a significantly excessive cross 
section at y/s = 3.773 GeV [TO] . The other method, al¬ 


lowing to consider the complicated interference effect, is 
to perform a scan around the 0(3770) resonance. Using 
this methd, the BESIII collaboration reports that the line 
shape of the cross section shows a deficit in the vicinity 
of the 0(3770) for the final states pp and ppn° Tlj T5] , 
Furthermore, there is a two-solution ambiguity for the 
branching fraction of 0(3770) -A pp/ppn 0 , which cannot 
be solved from the scan experiment. Recently, an evi¬ 
dence of 0(3770) —)• K + K~ was also found by studying 
the cross section of e + e“ —> K + K~ above 2.6 GeV [16j . 

We focus on the decays of 0(3770) going to baryon 
anti-baryon pairs (BB). Here B = p, A, X + , E°, and 
S°. All available experiments measuring the cross section 
of e + e _ —> BB at the c.m. energy from 3 GeV to 3.9 GeV 
are combined. In Sec. II, we will present the born cross 
section formulas of e + e _ —»• BB and introduce the model 
to relate all the baryon octet states. In Sec. Ill, we 
will review the available experiments and describe the fit 
strategy. The results will be shown and discussed in Sec. 
IV. A short summary will be given in Sec. V. 


II. CROSS SECTION FORMULAE OF e+e” BB 
AND DESCRIPTION OF THE MODEL 


The Born cross section of the QED process e + e _ —> 
7 * —> BB at the center-of-mass energy y/s can be written 
as 


vqed(s) = 


Attq. 2 Pb 

3s 


|G^(s)| 2 + ^|G|(s)| 2 


(1) 
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where Mb is the nominal baryon mass, /3 b = 

a/T - 4M%/s, a is the fine-structure constant, and 
and G B are the magnetic and electric form factors mi 
respectively. 

The resonance production cross section of e + e _ —>■ 
■0(3770) —> BB is written as 


i27rr e r b 

aRES = (s - Mo ) 2 + m 2 t 2 


( 2 ) 


where M 0 = 3773.15 MeV/c 2 and T 0 = 27.2 MeV [E] are 
the nominal mass and total width of 0(3770), T e (T0 is 
the partial width of 03770) —» e + e _ {BB). Tb can be 
written as 


the QED continuum amplitude. In many analyses (for 
example Ref. musi), this phase difference is assumed 
to be 0, namely, 0' = 0. We will find that the effect of the 
nonzero 0' is also negligible in the case of ■0(3770) —> BB. 

Therefore, the total cross section considering the in¬ 
terference between the processes e + e _ —>■ 7 * —>■ BB and 
e + e“ —>■ 0(3770) —► BB is constructed as 


<Jtot = 


Ana 2 /3 b 
3s 



2 M% C B 
s s 2 ln 2 (s/ A 2 ) 


r y/m/^y/T^E F B 
Vs s-M 2 + iM 0 T 0 \F B I 


( 6 ) 



where F B r and \F B \ are the form factors. 

The form factor ratio \Ge/Gm\ is 1 at the baryon 
pair threshold, but may have small deviations above the 
threshold. The predicted behavior is model-dependent 
(see for example Ref. mmi Experimentally, the form 
factor ratio is measured to be consistent with 1 within 
the uncertainties for the proton EB H23 in the region 
2.2 < 0s < 3.1 GeV and for the baryon A [21] in the 
mass region from the threshold to 2.8 GeV. However, the 
measurement of the neutron form factor from the thresh¬ 
old up to 2.44 GeV [21] indicates \Ge\ =0. Throughout 
this paper, we assume that \G B \ = \Gfj\ = \G B \ and 
\F B \ = \F^\ = \F B \ for all BB final states. The effect of 
this assumption will be considered. The nucleon electro¬ 
magnetic form factors in the timelike region have been 
extensively reviewed in Ref. [25]. Here, the form factors 
G b and F b take the following forms 


r^B 

G b {s ) =- % - 

V ' s 2 ln 2 (s/A 2 ) 


(4) 


To relate the form factors for all baryon octets, the 
SU(3) flavor symmetry is imposed. We also consider the 
SU(3) breaking effect due to the electromagnetic interac¬ 
tion and the quark mass difference of m s —m u /d- For con¬ 
venience, we introduce the matrix notations. The SU(3) 
octet baryons and anti-baryons are described by the ma¬ 
trices B and B respectively. 


B = 


/£ 0 /\/2 + A/06 
£- 


V 


£+ 

-E%/2 + A/V 6 

cto 


P 
n 

-2A/06) 
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B = 


/£%/2 + A/06 
£+ 


V p 


y- 

-£°/v / 2 + A/VQ 
n 



( 8 ) 

The SU(3) invariant effective lagrangian for the decay 
0(3770) —>■ BB can be written as 


Bff =gTr{ BB) 

+dTr{{ B, B}S e ) + jTr([B, B]S e ) 

+d'Tr{{ B, B}S m ) + f'Tr{[ B, B]S m ) , (9) 


from a calculation in Ref. [26] and 

F b = ]J^-^G B (M 2 )e^' + . (5) 

Here A = 0.3 GeV is the QCD scale parameter, C B and 
A B are the free parameters. In Eq. [ 5 J the first term 
represents the electromagnetic interaction amplitude of 
the 0(3770) and the second term represents the OZI- 
suppressed strong decay amplitude of the 0(3770). Two 
phase angles 0' and 0 are introduced relative to the QED 
process. 0' represents the phase difference between the 
electromagnetic amplitude of the 0(3770) resonance and 


where g,d,f,d', f are the coupling constants, “Tr” rep¬ 
resents the trace of a matrix, “[a,b]” and “{a,b}” denote 
the commutator and the anticommutator of the two ele¬ 
ments a and b respectively, and the matrices S e and S m 
are defined as 



In the right-hand side of Eq. [9] the first line represents the 
OZI-suppressed strong amplitude, the second line repre¬ 
sents the one-photon electromagnetic amplitude, and the 
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third line represents the SU(3)-breaking contribution due 
to the quark mass difference (more details about the ef¬ 
fective lagrangian can be found in Ref. [27H29] ). 

From Eq. [9j we can derive the following relations for 
the form factors G B and F B (or equivalently C B and 
A B ). 

C p = Ci + C 2 
C A = -Ci 
C s+ = Ci + C 2 
C E ° = Cl 
C H " = Cl - c 2 
C H ° = —2Ci 

and 


former two collaborations measure the cross section of 




■ pp using Eq. 11 


N° bs 


LiCi{ 1 + 5i) 


( 11 ) 


where i denotes the energy point, N° bs is the observed 
number of signal events, Li is the luminosity, ti is the 
efficiency, and (1 + Si) is the radiative correction fac¬ 
tor [33HSS]- The BABAR collaboration utilizes the ini¬ 
tial state radiation (ISR) technique [37]. The process is 
e + e _ —► 7 pp, where the photon can be required to be 
detected <22j or undetected [T!TI . The cross section of 
e + e _ — > pp at the c.m. energy of the pp invariant mass 


M pP can be extracted according to Eq. 12 


er(Mpp) = 


(dN/dMj 


pp)corr 


(. dL/dM p p)e(M p p)R 


( 12 ) 


A* = A 0 - Ai + A 2 
A a = A 0 - 2Ai 
A s+ =A 0 + 2Ai 

A^ = Aq + 2Ai 

A“ = Aq — A 1 — A 2 
A“ =■ Aq — A\ — A 2 

Here the free parameters Ci, C 2 , Aq, A\, A 2 are real num¬ 
bers in practice. 


III. EXPERIMENTAL REVIEW AND THE 
FITTING STRATEGY 

A. Experimental review 

We starts with the reaction e + e _ —> pp for which 
the most data sets have been accumulated. The BE- 
SIII collaboration has performed a scan from 3.65 GeV 
to 3.90 GeV and a deficit is found in the vicinity of the 
•0(3770) |14j. Considering the interference between the 
QED process and the 03770) resonant production, two 
solutions are found for the partial width of •0(3770) -A pp 
with equal goodness of fit. But Ref. [14] has not reported 
the statistical significance of the solutions. To solve this 
two-solution ambiguity, more experimental information is 
needed. The results from the studies of the proton form 
factors from the CLEO [HIED, the BES/BESIII [121 [32] 
and the BABAR [22) S3] collaborations can be used. The 


where ( dN/dM p p) corr is the mass spectrum corrected for 
the mass resolution effect, dL/dM pp is the ISR differ¬ 
ential luminosity [37], e(M pp ) is the detection efficiency, 
and R is the radiative correction factor. 

For the final states A A, E + E“, E°E°, S~S + and S°S°, 
the CLEO and BES/BESIII collaborations [TTMT3] have 
measured the cross sections at the peak of the 0(3770) 
resonance. Neglecting the interference effect with the 
QED process e + e _ — > 7 * — > BB, there is no significant 
excess compared to the cross section at an energy point 
far from any charmonium resonance. The BABAR col¬ 
laboration also studied e + e _ -a AA/E°E°/AE° using the 
ISR technique and provided the upper limit of the cross 
section at the 90% confidence level (CL), which will be 
used as a cross-check for our results. All the data sets 
used in the following fit are summarized in Tabic [I] and 
Table [TT] The denotations for the pp final state in the first 
column of Table |T] will be used consistently throughout 
this paper. 

B. The fitting strategy 

To combine the results from various experiments, we 
should consider the statistical uncertainties and the sys¬ 
tematical uncertainties correctly. For the number of sig¬ 
nal events N obs , it is either obtained by simply neglecting 
the background and counting the number of events or ex¬ 
tracted by subtracting the background events from the 
total number of events. Either way leads to a systemati¬ 
cal uncertainty. At all energy points in an experiement, 
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TABLE I. Review of the experiments of measuring the cross 
section of e + e _ —> pp from 3 GeV to 3.9 GeV. 


Denotation 
BES EH 

ESIII PH 


ip" scan 1141 


Process M pP (GeV) 
e + e“ —» pp 3.0, 3.07 

3.05, 3.06, 3.08, 3.40, 3.50 

e^e —> pp 

3.5507, 3.6002, 3.671 

3.65, 3.748, 3.752, 3755 

, _ _ 3760, 3.766, 3.772, 3.773 

e~”e —¥ pp 

3.778, 3.784, 3.791, 3.798 
3.805, 3.810, 3.819, 3.900 


BaBar (LA) [52] e + e — > 7 pp 3.0-4.0 
BaBar (SA) [33] e + e“ —¥ 7 pp 3.0-4.0 
CLEO EOIEI] e + e~^pp 3.671,3.772 


TABLE II. Review of the experiments of measuring the cross 
section of e + e“ —> BB from 3 GeV to 3.9 GeV. 


Denotation Process 


M BI3 (GeV) 


CLEO [TT] e + e —>■ A A —> ppn + n 


3.671, 3.772 


BESIII £3] 


E + E~ —> ppA-y 
E°S° —¥ ppn + n~ 2 7 
E _ E + —¥ pp2n + 2ir~ 
S°E 0 —» ppn + Tv~A'y 


3.65, 3.773 


7(AA) —7(pp7T + 7r _ ) 3.2-3.6 

■ 7(E°E°) —¥ 7(pp7r + 7r _ 27) 3.2-3.6 

■ 7(AE°) -7- 7(pp7r + 7r 7) 2.9-3.8 


the luminosities are measured using the same method, 
the signal events are selected using the same set of con¬ 
ditions, and the radiative correction factors are obtained 
in the same way. Thus the systematical uncertainties 
related to them are independent upon the energy point 
and will be considered by introducing a free normaliza¬ 
tion factor for each experiment. 


We starts with the case of proton. A \ 2 is constructed 
in Eq. 13 for each experiment except for the “ip" scan” 
experiement nu 


xl(p) = 


(1 -/a) 2 


ind. 


E 


(A° 6s - f a A,J 2 


(A N? b °)l t . + (Ae/e/ Q A ia )2 


(13) 

where a denotes the experiment, i a denotes the i-th en¬ 
ergy point for the experiment a, N obs is the observed 
number of signal events, A is the expected number of 
signal events and defined as A = crLe(l + S) or crLeR 
as indicated in Eq. 11 and Eq. 12 ( AN obs ) to t. is the 


quadratic sum of the statistical uncertainty of N obs and 
the systematical uncertainty due to the background sub¬ 
traction or neglecting the background events. Ae is the 


statistical uncertainty of the efficiency e determined from 
a limited MC sample. £ 2 nd is the quadratic sum of the 
systematical uncertainties which are independent upon 
the energy point. It includes the systematical uncertain¬ 
ties due to the consistent selection criteria at all energy 
points, the trigger efficiency, the reconstruction efficiency 
of charged tracks, the efficiency corrections as used in the 
BABAR measurements [22] S3] , the measurement of the 
luminosities, and the radiative correction factors. To con¬ 
sider these systematical uncertainties independent upon 
the energy point, the free normalization factor f a is in¬ 
troduced for each experiment. Here two things should 
be noted. One is that we do not consider the correla¬ 
tion of various selection conditions. The other is that we 
assume the form factors satisfy \Ge\ = \Gm\ and thus 
we do not consider the efficiency uncertainty due to this 
assumption (typically, the efficiency with |Ge| = \Gm\ is 
5% — 10% different from that with |Ge| = 0 ■’55]). 

For the “ip" scan” experiment in which N obs is found to 
be 0 at some energy points and the background contam¬ 
ination is only 0.6%, it is better to construct the likeli¬ 
hood function assuming that the number of signal events 
at each energy point abides by the poisson distribution 
as shown in Eq. [l4j 

X 2 p (p) = ( E 7 -2 ^ In P(N° b ° |/A0—In P(N° bs |lV° bs ), 

^>ind. ,• 

(14) 

where P(N\i/) is the probability of observing N events 
with the expectation value v in the poisson distribution, 
namely, P(N\u) = v N e~ v /AH, and / is the free normal¬ 
ization factor. 

For other baryon octets, the cross section of 
e + e _ —> AA at the peak of the ip( 3770) reported 
by the CLEO collaboration El and that of e + e —> 
E+E-/E 0 E 0 /S-S+/S°S 0 reported by the BESIII col¬ 
laboration m are used. As shown in the second col¬ 
umn of Table [TT] these processes share some final par¬ 
ticles such as protons, pions and photons. The related 
systematical uncertainties due to the reconstruction of 
proton and pion tracks, the second vertex fit, the parti¬ 
cle identification, and the detection of the photons are 
shared. However, Ref. [T3l did not report the individual 
systematical uncertainties. It is impossible to treat them 
correctly. Fortunately, the limited knowledge of the an¬ 
gular distribution contributes the dominant systematical 
uncertainty of 9.2% — 10.9%, which depends upon the 
baryon pairs and should be considered individually. The 
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X 2 is then constructed as follows. 


X\B) = Y, 

B= A,- 


(- N° B bs - A B? 

(AN° bs )l t . ± (AAb)L. 


(15) 


where B denotes the baryon, A is the expected number 
of signal events and defined as A = aLe(l + 5) x Bf 
with Bf being the product of the branching fractions of 
the intermediate-state decays, and (AA)tot. is the total 
uncertainty of the expected number of signal events. 

To combine all experiments, the full optimization 
quantity is defined as x%u = E a xlip) +Xp(p) + X 2 (B). 
Ao, 1 , 2 , 6 * 1 , 2 ) fa, ft and 0 are the free parameters. 


IV. FIT RESULTS AND DISCUSSIONS 

A. Fit to the cross section of e + e _ —> pp 


At first, we try the fit in the case of pp. The free pa¬ 
rameters are f a , A p , C p and </>. Here, ft is fixed to be 0 
for two reasons. One is that we can directly compare the 
result and that from Ref. M- The other is that float¬ 
ing ft leads to negligible difference. Two solutions are 
found with the same goodness of fit x 2 / n dof = 25.9/29, 
where ndof is the number of degree of freedom. The 
branching fraction of 0(3770) —> pp is found to be ei¬ 
ther (6.81™) x 10 ~~ 6 or (2.5 ± 0.1) x 10 -4 . If the pro¬ 
cess of 0(3770) -A pp is not included, the fit gives 
X ' 2 /ndof = 91.5/31, which means that the statistical 
significance of both solutions is larger than 5 standard 
deviations. Our results, summarized in Table |III| and 
Table IV are consistent with those in Ref. [14j . But 
Ref. Hi] does not report the statistical significance of the 
solutions. Figure [I] shows the cross sections of e + e“ —► pp 
from various experiments and the fit. There is an obvious 
deficit in the vicinity of the ■0(3770). 




vs (GeV) 


FIG. 1. The cross section of e + e - —> pp and the fit. (b) is the 
zoom of (a) around the 0(3770) resonance. Dots with error 
bar are the data. The black curve is the theoretic prediction 
from the fit. 


TABLE IV. Normalization factors from fitting to the cross 
sections of e + e _ —> pp 


Normalization factor Value 


fcLEO 

1.00 ±0.02 

fBabarS A 

1.01 ± 0.04 

fBabarLA 

0.99 ±0.05 

fBES 

0.99 ±0.05 

fBESIII 

1.00 ±0.03 

fft}"scan 

1.02 ±0.05 


B. Fit to the cross sections of e + e —> BB 


TABLE III. Two solutions from fitting to the cross sections 
of e + e“ —> pp 


B(0(377O) -A pp) 

A p 

C p (GeV 4 ) 0 

(6.8 tH) x 10- 6 

1 4+0.6 
i * 4 -0.3 

61.7 ±1.5 (-109° ±28.3°) 

(2.5 ±0.1) x 10~ 4 

8.4 ± 0.2 

61.7 ±1.5 (-94.4° ± 5.3°) 


Including all experiments, the fit results are summa¬ 
rized in Table |V| and shown in Fig. [2| The goodness of 
fit is x' 2 /ndof = 43.7/31. From Fig. I we find that the 
line shape shows a dip structure around the 0(3770) res¬ 
onance for the final states pp and £ + E - and a bump 
structure for the final states AA, E°E°, and S°S°. 




































6 






FIG. 2. The line shape of e + e“ —► BB and the fit. (b) and (d) are the zoom of (a) and (c) around the ?/>(3770) resonance, 
respectively. Dots with error bar are data. The curves are the theoretic prediction from the fit. In (a) and (b), the color-baryon 
correspondence relation is black-p, red-A, green-E + , blue-E°, yellow-E _ and pink-E°. The measurements of <r(e + e~ —> pp) are 
drawn in (c) and (d) while omitted in (a) and (b) for convenience. 


TABLE V. Fit results for e + e —» BB 


Parameter 

Value 

A 0 

2.91 ±0.22 

Ai 

0.37 ±0.15 

A2 

0.07 ±0.30 

Gi (GeV 4 ) 

19.6 ±2.3 

C 2 (GeV 4 ) 

41.7 ± 2.6 

4>' 

-54° ± 230° 

<t> 

-137.5° ± 2.7° 


In Table |VJ we find that | ^4i |, |^4 2 1 << |^4o I > which 
means that the SU(3) breaking effect is small. In addi¬ 
tion, the upper limits of a(e + e~ —> AA/E 0 E°/AE 0 ) at 
3.2 — 3.6 GeV from the BABAR measurement [23] are 
consistent with the predicted cross section from the fit 
result. Using the parameters from the fit, the branching 
fractions of ^>(3770) -± BB are calculated according to 
Eq. [ 3 ] and listed in Table VI All B(ip( 3770) — >■ BB) s are 


of the order of 10 5 . 


TABLE VI. Branching fraction of f/>(3770) —> BB. The first 
uncertainty is from the global fit and the second uncertainty 
is due to the assumption that the electric and magnetic form 
factors are equal. 

Baryon B(ip( 3770) BB) (xlCT 5 ) 


pp 

2.4 

± 

°o 

0 

± 

CO 

0 

AA 

1.7 

± 

0.6 

± 

0.1 

E+E“ 

4.5 

± 

0.9 

± 

0.1 

E°E° 

4.5 

± 

0.9 

± 

0.1 


2.0 

± 

0.7 

± 

0.1 

0 

mi 

3 

to 

2.0 

± 

0.7 

± 

0.1 
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C. Some discussions 


1 . In the analysis above, the relaitons \Ge\ = \Gm\ 
and \Fe\ = \Fm\ are assumed. Fits are repeated 
assuming \Ge\ = 0 and \Fe\ = 0 instead, which 
is indicated from the measurement of the neutron 
form factors [24] . The branching fraction difference 
is taken as the systematic uncertainty (the second 
uncertainty term in Table VI). 


2. We find that the two-solution ambiguity of 
£>(0(770) —>■ pp) reported in Ref. P2Q is fixed with 
including the measurements about other baryon 
pairs. This can be clearly shown by comparing the 
y 2 curves as a function of the parameter A p using 
only the cross sections of e + e _ —> pp and using the 
cross sections of e + W —>• BB (B = p, A, E + , E°, S~ 
and E°). The reduced y 2 curves are illustrated in 
Fig-! The reduced y 2 is defined as the difference 
of the y 2 from the fit with A p fixed and that from 
the best fit. The blue curve in Fig. [3] indicates that 
the smaller solution in Ref. m gives a better fit 
including all measurements of cr(e + e _ —>• BB). In 
the case of assuming \Ge\ = 0 and \F E \ = 0, this 
conclusion does not change. 



FIG. 3. The curves of the reduced y 2 as a function of A p 
using only the cross sections of e + e“ — > pp (red curve) and 
using the cross sections of e + e“ —> BB (blue curve). The 
reduced y 2 is defined as the difference of the y 2 from the fit 
with A p fixed and that from the best fit. The arrows denote 
the best fits. The solid and dashed horizon lines denote the 
lcr and 2cr regions respectively. 


within the uncertainty. Fixing <p' = 0 produces 
negligible effect. The relative phase of the OZI- 
suppressed strong decay amplitude of the 0(3770) 
and the QED amplitude is found to be —137.5° ± 
2.7°. Many phenominological analyses [251 HTSffTo] 
have been performed for various final states in the 
hadronic decays oi J/ip and 0(3686). It is revealed 
that the relative phase (p is close to —90°. For 
J/ip/ip{ 3686) —> BB, two possible phase values are 
found using a similar model in Ref. [29]. If the 
relative phase is assumed to be universal whatever 
the final state is, the large negative values are fa¬ 
vored and close to —90° for J/ip and ^(3686). For 
the decay mode 0(3770) —> pp, the calculation of 
Ref. [45] shows that the dominant contribution is 
the OZI-suppressed amplitude and (p = —113°. If 
the contribution of the OZI-allowed DD state as 
an intermediate state (which is firstly introduced 
in Ref. [47]) is included, the phase angle (p becomes 
—99°. However, our study shows the phase is far 
from —90°, which indicates that there may be ad¬ 
ditional mechanism contributing to the baryon-pair 
decays of 0(3770). 


TABLE VII. The relative phase between the amplitude of the 
strong interaction and that of the electromagnetic interaction 
in the decays of ip —> BB. 


Charmonium 

0 ' 


0 

J/iP 

fixed at 0 

-85.9° 

± 1.7° or +90.8° ± 1.6° [29] 

0(3686) 

fixed at 0 

1 

co 

00 

o 

± 25° or +134° ± 25° [29] 

0(3770) 

-54° ± 230° 


-137.5° ±2.7° 


4. Assuming the probability of a charmonium state 'ip 
decaying to light hadrons is proportional to the ab¬ 
solute square of the value of the charmonium wave 
function at the origin, we can relate the branch¬ 
ing fractions of J /ip/ip (3686) /ip (3770) —> BB. We 
define a dimensionless quantity, k^, in Eq. 


b F(0 t BB)/f3* B(ip BB)/P. B 

* “ r M) Be(ip) ’ ( j 


16 


3. The relative phase between the electromagnetic 
amplitude of the -0(3770) and the QED ampli¬ 
tude (t p' = —54° ± 230°) is consistent with 0 


where 1 - 4M 2 /M 2 and T e (ip) ( B e {ip )) 

is the partial width (branching fraction) of ip —> 
e + e~. Under the assumption above, it is expected 




















that Kji^ = k. 


B 

0(3686) 


= K. 


B 

ip( 3770) 


and that 


0(0(3770) -A SB) = 


K 


B 

03770) 




B sc .(ip(3770) -A BB) = 0(0 


0sc. (0(3770) -A BB) (17) 

^ BS) BeM 3770 )) 


W) 


/3.f 


where ip = J/ip or -0(3686), and i5 sc . (0(3770) —► 
BB) is from a scaling of 0(0 —► BB) under the 
assumption k^( 377 o) = 

lists the fi^s for ip = «//-0, -0(3686) 


VIII 


Table 


and 0(3770). 


We find that Kj/^ — 


K 03686) < 


0.1k^( 37 7o), which means that 0(0(3770) -A BB) 
is at least one order of magnitude larger than that 
scaled from B(J/ip/ip (3686) -A SB), as shown in 

Eq.irn 


TABLE VIII. k values for 0 —>■ BB. For J/ip and 0(3686), 
the experimental measurements included in the PDG [TH] are 
used. 


Baryon pair Kj/j, 



*±(3686) 


*±(3770) 

PP 

0.045 

± 

0.001 

0.041 

± 

0.002 

2.9 

± 

0.9 

AA 

0.039 

± 

0.004 

0.045 

± 

0.008 

2.2 

± 

0.8 

E+E~ 

0.039 

± 

0.006 

0.043 

± 

0.013 

6.1 

± 

1.3 

E°E° 

0.034 

± 

0.002 

0.037 

± 

0.007 

6.1 

± 

1.3 


0.027 

± 

0.004 

0.033 

± 

0.011 

2.9 

± 

1.0 

[I] 

o 

[I]l 

O 

0.024 

± 

0.003 

0.051 

± 

0.016 

3.0 

± 

1.0 


TABLE IX. Partial width of the observed non-DD decay 
modes of the 0(3770). r mens . (Y pr ed.) represents the mea¬ 
sured (predicted) partial width. 


Decay Mode Y meas . (keV) 

r pred. (keV) 

1 

+ 

52.5 ±7.6 

20-110 gg 

o 

o 

21.8 ±8.2 


J/07? 

24 ± 11 


7XcO 

198.6 ± 24.5 

199-225 gHEE], 524 gl] 

7Xd 

73.4 ± 13.6 

59-77 [HH52] 

e + e“ 

0.262 ±0.018 


H 

8.4 ± 1.9 



However, for exclusive light hadron decay modes, 
it is difficult to have an accurate theoretical predic¬ 
tion. In this work, the combined branching fraction 
of ip( 3770) —► BB is of the order of 10 -4 by sum¬ 
ming up the numbers in Tabic |VI| Though it is 
much smaller than the non-BB branching fraction 
of the order of about 10% measured by the BES col¬ 
laboration PEI, the baryon pairs only account for 
a small fraction of the light hadron decay modes. 
Furthermore, this work shows B(ip( 3770) —► BB) 
is at least one order of magnitude larger than that 
scaled from B(J/ip/ip( 3686) —> BB) as discussed 
above. This indicates that the mechanism in the 
light-hadron decays of the -0(3770) is different from 
that in the case of the J/ip or 0(3686). 


5. In view of last point, it is necessary to have a small 
review about the non-DB decay modes which have 

summa- 


been observed experimentally. Table IX 


rizes the measured partial width of these non -DD 
decay modes and the theoretical predictions. The 
measured partial width is calculated by multiply¬ 
ing the full decay width by the corresponding mea¬ 
sured branching fraction [T 8 ]. From Table IX 


we 


find that the potential models proposed in Ref. [48] 
and Ref. [1§1 - (KT] can explain well the rate of the de¬ 
cay modes with the charmonium final state. These 
models assume that 0(3770) and 0(3686) are the 
mixture of the 2 S and IB states of the cc system, 
as shown in Eq. 18 with the mixing angle 9 ~ —10°. 


0(3686) = 02 s cos 9 ± ipuj sin 8 

0(3770) = —ip 2 S sin 9 ± ip 1D cos 9 (18) 


V. SUMMARY 


Focusing on one type of non-BB decays, 0(3770) 
into baryon anti-baryon pair, all available experiments 
of measuring the cross section of e + e” —► BB at c.m. 
energy ranging from 3.0 GeV to 3.9 GeV are collected. 
A model based on the SU(3) flavor symmetry is built to 
relate the baryon octets. The SU(3) breaking effects due 
to the electromagnetic interaction and the quark mass 
difference are also considered. A global fit based on this 
model is performed. The two-solution ambiguity about 
0(0(3770) -A pp) reported in Ref. [2] is fixed. We find 
that the statistical significance of the presence of the 
process e + e~ -A 0(3770) -A pp is much larger than 5 
standard deviations, which is not reported in Ref. M- 
0(0(3770) -A BB) is determined to be (2.4±0.8 ±0.3) x 
10~ 5 , (1.7 ± 0.6 ± 0.1) x 10 -5 , (4.5 ± 0.9 ± 0.1) x 10" 5 , 
(4.5 ± 0.9 ± 0.1) x 10~ 5 , (2.0 ± 0.7 ± 0.1) x 10" 5 , and 
(2.0 ± 0.7 ± 0.1) x 10" 5 for B = p, A, S+, £°, 5" and 3°, 
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respectively. They are at least one order of magnitude 
larger than a simple scaling of $( J/r/>/-0(3686) —► BB). 
Furthermore, the relative phase between the strong am- 
pltitude and the electromagnetic amplitude is found be 
to far from —90°, which are favored in the hadronic de¬ 
cays of J/ip and -0(3686). The two evidences above may 
shed light on the puzzle about the non -DD branching 
fraction of ^(3770). 
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